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We consider active suspensions in the isotropic phase subjected to a shear flow. Using a set 
£\J \ of extended hydrodynamic equations we derive a variety of analytical expressions for rheological 

quantities such as shear viscosity and normal stress differences. In agreement to full-blown numerical 
calculations and experiments we find a shear thickening or -thinning behaviour depending on whether 
the particles are contractile or extensile. Moreover, our analytical approach predicts that the normal 
I , stress differences can change their sign in contrast to passive suspensions. 
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<£J I. INTRODUCTION 

o 

Within the last years the physics of active materials has become a focus of growing interest. In an active system, 
each particle (constituent) " constantly absorbs energy from its surroundings or from an internal fuel tank (such as 
adenosinetriphosphate) and dissipates it in a process of carrying out translational or rotational motion" [8]. The 
resulting, continuous energy flux drives an active system out of equilibrium even in a steady state. This is in sharp 
contrast to the physics of "passive" soft- matter systems where non-equilibrium behavior typically occurs only in the 
£h ' presence of an external driving force (such as shear) . 

Active systems often occur in biological contexts, examples being acto-myosin filaments, cytoskeletal gels interacting 
i with motors, bacterial colonies, and swarming fishes [H-0|. Moreover, there are also non-biological representatives such 

as vibrated granular rods 0, H| or remotely powered, self-propelled particles @. 

The individual motion of an active particle through the suspension gives rise to flow fields which, on long time 
scales, can be modeled by force dipoles n 0, H|. The combination of this directed motion and interactions between 
active particles (e.g., excluded volume interactions or hydrodynamic interactions) generates rich collective behaviour 
with some similarities to that of ordinary "passive" liquid crystals. However, active systems ("living liquid crystals" 
@) have fascinating additional features such as an instabilit y o f homogeneous, nematic states or the occurrence 
of spontaneous flow states in films above a critical thickness jlll - llq ] . 

Modeling these phenomena on a microscopic level is still a challenge Q . As an alternative, the collective behavior is 
often studied via a coarse-grained continuum approach familiar from liquid-crystal theory. Within this framework the 
relevant dynamic variables are the tensorial orientational order parameter = (n^nj) [ljj (where describes 
the particle orientation and tt7 indicates the symmetric traceless part of a tensor) , the velocity field, the stress tensor 
<7 MI ,, the density (if the latter is allowed to vary), and optionally, the polarization of the system (ll M. flil \l$L I2H I22I] . 
The latter becomes particularly important in systems of so-called "movers" such as fishes and bacteria. The resulting 
hydrodynamic equations are then supplemented by additional terms which take into account the active nature of 
the system (see, e.g., Refs. [1, [TH, [I3])- As a minimal ansatz it is often assumed that the individual force dipoles 
generate a contribution to the stress tensor which is linear in Q^. The prefactor can be either positive or negative 
describing, respectively, "contractile" systems (such as suspensions of alga Chlamydomonas Reinhardii, cytoskeletal 
rods, actomyosin) or "extensile" particles (e.g., E. Coli bacteria). 

In the present work we use a hydrodynamic approach to explore the non-linear response of an active material to an 
external shear flow. We focus on a non-polar material composed of " shakers" (such as melanocytes or symmetric rods) 
and neglect any variations in the number density. The passive contributions to the equations of motion are modeled in 
the spirit of our earlier work on passive liquid-crystalline systems (see, e.g., [23|-[25j]). Activity is then incorporated by 
the above-mentioned stress term plus the lowest-order contribution to the order parameter itself. Previous theoretical 
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investigations on the rheology of active liquid-crystals have already indicated that their flow behavior is more complex 
than that of the corresponding passive system. On top of the above-mentioned spontaneous flow one observes, e.g., 
activity-induced shear band formation [3, E3] , hydrodynamic instabilities 0, [l(| , fluid mixing [2(| , pattern formation 
in steady states [13, HI] and cell-sorting phenomena [29| . 

Contrary to many of these studies, which often involve quite complex numerical investigations, we aim to derive 
analytical results. The focus is on experimental accessible quantities such as the shear viscosity and the normal 
stress differences. Our motivation to approach these quantities from an analytical side is to establish simple "rules" 
for the role of, e.g., the prefactor of the activity- induced pressure term, on the non- Newtonian viscosity. We hope 
that the physical insight gained by such simple rules may be particularly interesting for the interpretation of future 
experiments. Moreover, the results can also serve as benchmark for more complex numerical studies. Of course, 
the derivation of explicit expressions is impossible for the full, nonlinear set of hydrodynamic equations. The main 
assumption in our work is that the system is close to the passive isotropic state. From a physical point of view this 
means that we are essentially restricted to the investigations of states which are isotropic in passive equilibrium in the 
absence of external flow. Despite this assumption, our theory shows a variety of non-linear flow phenomena such as 
shear thickening and shear thinning in dependence of the activity parameter. Moreover, the theory yields novel results 
for the normal stress differences, which are a characteristic rheological feature of any soft-matter system beyond the 
linear-response (Newtonian) regime. 

The remainder of this paper is organized as follows. In section [II] we present our set of hydrodynamic equations and 
motivate (following the arguments in earlier studies) the extra terms due to activity. Section IIIII deals with the shear 
viscosity (in a slab geometry), focussing on a spatially homogeneous system. Here we present analytical expressions 
and discuss our results in the light of earlier, fully numerical studies. Section IIIIBI is devoted to the normal stress 
differences on the same level of approximation. 



II. MODEL EQUATIONS 

In the continuum limit the orientational dynamics of active particle suspension is governed by a coarse-grained 
equation of motion for the second-rank order parameter Q^ v Q . The starting point are the well-established relaxation 
equations for passive liquid crystals pBI ] , which may then be supplemented by terms reflecting the activity. Including 
these effects to lowest order [H, , the resulting equation for the time derivative of Q^ v is given by 

(d t + waVa)Q M! / = e^xn^xQitu + "^kT^xQxu - V^Tq 1 TQ P T liV 

-Tq 1 ^ + — AQ M , + AQ^, (1) 
T Q 

where we use the Einstein summation convention and the notation for the velocity. The first three terms on the 

right-side of Eq. (fTJ) describe the impact of the flow on the order parameter. Specifically, V ^ = V^u^ is the strain 
rate, w^jthe vorticity, and tq p and tq are relaxation time coefficients, respectively. In the subsequent term in Eq. ((T|), 
= is the "molecular field", which is described by the derivative of the Landau-de Gennes (LG) free 

energy $ of a homogeneous, liquid-crystalline system, 

$ = \A(T)Q^ V Q^ - J^BQ^Q Xv Q^ + - A C{Q^Q^) 2 , (2) 

where A(T) — a(T—T c ) (with T being the temperature) changes its sign at the isotropic-nematic spinodal temperature, 
T c . Note that the definition of a free energy of an active system is, strictly speaking, impossible due to its non- 
equilibrium character. We thus consider $ as the system's free energy in the "passive phase". The next (gradient) 
term in Eq. (TTJ) is again borrowed from equilibrium physics; it results from the increase of free energy due to distortions 
in the one-constant approximation. The prefactor of this "elastic" contribution corresponds to the squared elasticity 
length, £q. So far, the terms in Eq. ([T]) correspond to those for passive systems (see, e.g., 25]). It is the last term in 
Eq. ([1} which corresponds to an ansatz for the effect of activity on the orientational order, with A being a coupling 
parameter. For A > (A < 0) the order parameter increases (decreases) in time; the case A > is therefore 
sometimes referred to as " self- aligning" . For a microscopic interpretation of the coefficient A, see (3f|. 

Some further interpretation arises when we consider this term, which is linear in Q M „, together with the lin- 
ear term resulting from the derivative of the squared term in the LG free energy. This yields a contribution 

^— Tq 1 j 4(T) + Xj Q M „ = — TQ~ 1 S$/SQ f _ lUl where <l is a pseudo free energy defined by 

® = \{A{T)-\r Q )Q llv Q llu . (3) 



3 



From this consideration we see that A influences the onset of nematic ordering upon cooling the system from high 
temperature (i.e., large values of A). Specifically, a positive value of A effectively increases the spinodal temperature 
and thus stabilizes nematic ordering with respect to the passive case. On the other hand, A < stabilizes the 
isotropic (i.e., orientationally disordered) state We also see from Eq. ([3]) that for A > Tq 1 A(T) the isotropic state 

is unstable (i.e., 8 2 §/(5Q) 2 < 0). Finally, Eq. ([3]) shows that the order of magnitude of A is given by the inverse of 
the relaxation time tq . 

The fluid velocity field obeys the generalized Navier-Stokes equations 

p(d t + «aVa)k p = ZrHsodxTxfj, + d x ax^,, (4) 

where p is the fluid density, and ?7i so is the first Newtonian viscosity of the passive system. Following earlier studies 
(see, e.g. [ToL Il5j|) we model the stress tensor o~xn as a sum of a passive part, <r p , and an active part, <7 a , i.e., 

The passive part is given by the corresponding liquid-crystal expression [23| 

= --fcBlW^P^- AQ M „) (6) 
m \tq 

- xGkQuxQxJ + V2k $ 'Q&AQjJ) , (7) 



which can be derived from the free energy of the system and the principles of irreversible thermodynamics. The 
presence of the active term, cr a , results from the fact that the self-propelled particles in an active system induce flow 
fields which are dipolar in character [H, [l(|. Here we use the simplest ansatz (for higher-order corrections see [3ll |) 
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-k B T(Q^. (8) 



The proportionality constant distinguishes between extensile (£ > 0) and contractile (£ < 0) active flow behavior 
H O E3| and is characterized by the strength of the elementary force dipoles. From Eqs. (jHJ) and (jlj we then see 
that an extensile system opposes the external flow (described by T,,^), whereas a contractile system enhances it. 

Finally, our analysis neglects any fluctuations of the concentration of active particles in the suspension. Clearly 
there are conditions, where the dynamics of the concentration field can strongly influences the other variables 16] . The 
goal of the subsequent analysis is to investigate the influence of the parameters A and C on experimentally accessible 
rheological properties of the active system. In doing this we focus on the passive, isotropic state rather than on the 
metastable regime at the isotropic-nematic transition which was studied numerically in Ref. |T7| . 



III. ANALYTICAL RESULTS 



The flow geometry is given by two infinitely extended plates in the x-z-plane, which are separated by a distance 
2h along the y direction. We apply a planar Couette flow by moving the upper plate by +v w and the lower by — v w , 
yielding a velocity profile u M = (v(y),0, 0). Under these conditions the strain-rate tensor simplifies for a linear flow 
profile to V ^ v — d y v x /2S fiy S I/x = j /28^ y 8 vx , where 7 = dv/dy is the shear rate. 

To rewrite the tensorial equations of motion Eqs. (flj and (|6]) in a more explicit way, we express the second-rank 
tensors Q^ v and a^, which have five independent components, in a tensor basis, e.g., — X)a=o QsXau [23 • The 
basis tensors are defined by 

T° — \ -e z e z ' T 1 — \ -<e x e x ~ e y e y ) T 2 — \/2e x e y ' 
± fif ■— y 2 /* f iu> '~ y 2 ^ P v ii v)i ytv '~ v ^n^v 

1% := V2<^, 1% := V2efX, (9) 

where, e.g., is the /^-component of the unit vector in z-direction. 

In the subsequent analysis we assume that the passive system is isotropic in the absence of external flow (i.e. 
Qfj,u = 0). In this case the order-parameter expansion in the equilibrium free energy can be truncated after the first 
(quadratic) term, yielding only the linear contribution (i.e., <^^ v sa AQ^) in the dynamic equation for Q^. We 
also assume, as a first approximation, that the system is spatially homogeneous, that is, AQ = 0. The resulting, 
homogeneous equation of Eq. ([1] ) for the alignment tensor decouples into equations for the symmetry-adapted 
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(a= 0, 1, 2) and symmetry-breaking (a= 3, 4) components. The symmetry-breaking components give no contributions 
to our rheological quantities. Therefore, we focus on the symmetry-adapted components, that is, 

Qn» = QoT° v + Q x Tl v + Q 2 T^ . (10) 

A similarly decoupling occurs for the stress tensor, yielding 

eV„ = ^oT° y + aiT^ u + a 2 T% u . (11) 

Inserting Eq. (ITU|) into Eq. ([T]) we obtain the following set of equations describing the orientational dynamics in a 
homogeneous, isotropic active system: 

TQ 1 + -^=kTQ q + Q 2 - -^r + A^Q 2 
V3 tq A 

Qi = TQ 2 + \^Q l 

Qo = Q 2 + a^q , (12) 

where T = tq^/A. Equations (fl"2"|) can be solved for the three symmetry-adapted components yielding, e.g., 

Q 2 (r) = -^p .£ r . (13) 

t q i — a* + (i — ^K 2 )r 2 (i — a*) 

Similar expressions result for Qq and Q\. In Eq. (| 13|) we have introduced the dimensionless active parameter A* = 
Atq/^4 related to the aligning term in Eq. ([1]). Interestingly, there is no dependence of Q 2 on the other parameter 
related to activity, £* = (/A, which appears in the pressure equation. We note that this is not a consequence of the 
linearization, but rather of our assumption of spatial homogeneity (no backflow). 

For liquid crystal polymers the shape parameter k depends on the coupling strength between orientation and fluids 
strain rate. In the kinetic approach n can be related to the relaxation times tq p and tq |36| . i.e. k = — \/15tq p /(7tq). 
Note, in the description of passive liquid crystal polymers the entropy production has to be positive and Onsager's 
symmetry relations apply. As a result r p , tq > 0, but tq p can have either sign depending on the particles shape. For 
rod-like particle suspensions the relaxation time tq p is negative and for disk-like particles is positive, respectively. For 
active suspensions Onsager's symmetry relations may be broken such that tq < is possible. However, for simplicity 
we only discuss tq > 0. 

Some other features of Eq. (|13|) can be seen immediately. In the limit T — > (no external shear flow) Q 2 vanishes 
(as do the other components), consistent with our assumption of a passive isotropic equilibrium system. A special 
situation occurs for A* — > 1, which corresponds to approaching the stability limit of the isotropic phase [see Eq. ©]. 
At the same time, the order parameter as determined from Eq. (QJ increases to infinity. This behavior contradicts 
the physical picture according to which the order parameter is restricted to a finite (saturation) value related to 
perfect alignment. Therefore, the limit A* — > 1 corresponds to a situation where the truncation of the LG free 
energy after the quadratic term cannot justified any more, i.e., higher order terms need to be included to guarantee 
the existence of a minimum. Later, we will see that this divergence gives rise to a divergence of the viscosity for 
both extensile and contractile suspensions. A divergence for contractile suspensions was predicted by Marchetti and 
Liverpool [1, HI, US HH an d confirmed in numerical simulations by Cates et al. pTj . 

For extensile suspensions the viscosity drops to zero [ItJ ■ As we will later see, a different behavior is predicted 
by our model where the viscosity diverges to minus infinity. This unphysical result effectively restricts our study to 
values of A < 1. 



TABLE I: Parameters 

TQp T p TQ K 

-0.1 0.1 0.2 0.37 



A. Apparent shear viscosity 



The rheological properties of non-Newtonian fluids can be captured experimentally by measuring the shear viscosity 
and the normal stress differences. In the geometry used here, the shear viscosity is determined by the xy— component 
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FIG. 1: (Color online) The shear stress vs. shear rate T for A* = 0.7. The remaining parameters are given in table U 
Left: Extensile suspensions show a plateau and negative stress close to zero shear rate. Right: Contractile suspensions are 
characterized by a yield stress. 



of the stress tensor, i.e. 5™ = v2<72. I n terms of the shear rate V the shear stress reads 



r Qp, 



n- 2 ;; ls(1 , - -k B TA ( ^(1 + -V :Q, 

m \ tq r Qp 



3 1 - A 



7Q2 



(14) 



where ry; so = (p/m)kBTT p (l — t q p / ( t Q t p)) is the Newtonian viscosity in a system with Q = (i.e., in an isotropic 
state or for systems of spherical particles) . 

The apparent shear viscosity 77 is determined by the ratio of shear stress versus the shear rate. For our analysis we 
scale 77 with the (passive) first Newtonian viscosity ijncw = p/mkBTT p , yielding 



= 1 



'Qp 
T Q T P 



C* 



1-A* + (1- |k 2 ) x 



TQpl-A* + (l-i K 2 )T rjL 



(15) 



For vanishing activity parameters A* = and £* = the apparent shear viscosity of a passive liquid crystal is recovered 
23]. 

In Fig. [TJ the flow curve of our model is shown for extensile and contractile suspensions. For extensile suspensions 
there exists a threshold Q = 0.55 above which the shear stress shows a plateau. A plateau in the constitutive curve 
is frequently used to explain shear banding flow instabilities in passive complex fluids that have been studied for a 
long time [19j . The difference here is that the plateau occurs around the zero shear point and the related instability 
might generate shear bands with velocities directed in the opposite direction. At the same time the viscosity vanishes 
reminiscent to a superfluidic state as discussed in fvA [T§j . 

In a second possible scenario, at a critical nonzero stress — 00 (ay) the shear rate jumps from a positive (negative) 
to a negative (positive) value forming a hysteresis [l8j . Hysteresis effects around a positive shear rate value can be 
observed for anisotropic complex fluids [2fJ. In contrast to the passive system for active suspensions the hysteresis 
enclose the zero shear rate point. This means that in a shear stress controlled experiment starting with a > <jq there is 
a threshold before the direction of the velocity can be flipped. A similar effect is not observed in passive suspensions. 
The threshold is related to the extra active stress due the active particles orientation and interaction to the solvent. 
Contractile suspensions do not show a plateau, but a yield stress which becomes more and more pronounced for higher 
values of C* • 

Passive liquid-crystalline systems typically display shear thinning, that is, the shear viscosity decreases with in- 
creasing r. This is indicated by the red dashed line in Fig. [5J The shear thinning results from the coupling of the flow 
gradient and the orientational degrees of freedom, yielding flow- alignment of the suspension. In our model all three 
components of the tensorial order parameter are effected. The eigenvectors and eigenvalues indicates the biaxiality 
of the flow aligned state. The arrows above Fig. [5] display the eigenvectors multiplied by its eigenvalues. A small 
angle is enclosed between the principal director (the eigenvector corresponding to the highest eigenvalue) and the 
flow direction. The lower part of Fig. [2] shows the viscosity for extensile and contractile suspensions. Obviously, the 
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FIG. 2: (Color online) Top: The components of the order parameter Qi for parameters given in[I] The arrows above the figure 
show the eigenvectors of Q. The length is given by the magnitude of the eigenvalues. 

Bottom: The scaled viscosity rf vs. shear rate V for £* = —0.6, 0.0, 0.6 and A* = 0.7. The remaining parameters given in table 
|U The red dashed line indicates shear thinning for the passive suspension. 



behavior if rj depends very strongly on £ at low shear rates, but becomes essentially independent at high shear rates. 
The intermediate regime is characterized by shear thinning or -thickening behaviour depending on the actual value of 

C- 

For passive liquid crystal polymers in the zero shear limit the viscosity is equal to the the Newtonian viscosity ?7i so . 
However, for active suspensions the zero shear viscosity also depends on active parameters, that is 

2 /A* + Tq C* \ 

*-i-^(-s^i-J- < 16) 

It is remarkable that in the active case the shear viscosity at T = strongly differs from the Newtonian case, indicating 
the non-equilibrium character of active suspensions. 

For contractile suspensions the shear thinning is enhanced and diminished for extensile. For extensile suspensions 
there is transition to shear thickening when the zero shear viscosity becomes smaller than the second Newtonian 



7 



viscosity 77^ = 1 — Tq /(tqt p ). For very high extensile strength the zero shear viscosity can vanish related to the 
plateau of the shear stress rate curve Fig. PR The effect of reduced viscosity was also discussed for a two-dimensional 
model of dilute bacterial suspensions in 138]. As mentioned before the high shear rate viscosity reaches the second 
Newtonian viscosity independent on the activity. 

In our model the shape parameter k influences the results only marginally. In the special case K = the shear 
viscosity simplifies to 

77* = 77* + 770 ~ 7,00 (17) 

where r r is a characteristic time defined as t^ 1 = A/tq — A* = 6D r — A*. Therefore, the activity parameter A* modifies 
the self-rotational diffusion constant D r . The apparent viscosity Eq. ( ITT]) can be related to a simple rheological shear 
thinning model, that is, the Cross model (46[: 

V = Voo+ 1 + {c . )m . (18) 

For our model we find C — r r and m = 2. 

The simplified formula (17) gives us the possibility to link experimentally accessible quantities to our theoretical 
parameters A* and £* . Specifically, we need from experiments the zero-shear viscosity, the high-shear viscosity 7700 , the 
rotational diffusion constant, and the relaxation time r r . The latter can be determined from a fit of the experimental 
shear viscosity to Eq. (16). The rotational diffusion constant D r has to be measured in an independent experiment, see 
e.g. [HI for rod-like fd- viruses. From these quantities, the activity parameter A follows via the relation A* = 6D r — t^ 1 . 
Next, the active force acting on the fluid modelled by (* is related to A* via 

^- C * = (r?r-1)A*-^, (19) 
T Qp 

where we defined the viscosity r/l — . Here the ratio of the relaxation times can be related to the active particles 

shape g3, i.e. ^ = -J^R- The coefficient R measures the non-sphericity of particles. For isolated ellipsoids with 

the semi-axes a — b, c and the aspect ratio q — c/b, one has R = (q 2 — l)/{q 2 + 1). 

There is a remarkable symmetry between the particles shape and response to the flow. The rheological properties 
do not distinguish between extensile disk-like particles (contractile, rod-like particle) and contractile rod-like particles 
(extensile, rod- like particle), respectively. 

B. Normal stress differences 

The appearance of normal stress differences indicate strong non-Newtonian behaviour and is related to surprising 
rheological effects, e.g., the Weissenberg effect or swelling jets [46| . In our theoretical description we can compute the 
normal stress differences from the components of the stress tensor via N± = a xx — a yy and N 2 = <J yy — a zz . In our 
notation 

#1=2(71, N 2 = — \/3(7o - (7i. (20) 
The analytical expression of the stress tensor components are given by 
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ID 



V3 V 3 ; (1-A*) 2 

From these relations it follows that 

/ 



1 + f2 M \ ) j (22! 



Ni_ = 2 jQ P A 



2 

^Ncw T Q T P 



A * + (1-| k2)t ^ 



2 



TQp 1- A* + (1-^)^ 



A* 



1 - A* 



(23) 
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Similarly N2 is given by 



N 2 



t Q t p V 



t Qp S 



1 - A* 



1-A* + (1 + |^) I ^ 



1-A* 

r 2 



r- 



(24) 



1-A* 



(1 



3 k ; X _ A . 



1-A* 



In passive, nematic liquid-crystal polymers subject to a shear flow the first normal stress difference displays a different 
sign depending on the shear rate. The resulting change of sign was first observed in experiments by Porter and Kiss 
[39l | (see also [4(J EH, EH for further experiments and numerical simulations) . 

On the other hand, isotropic liquid crystal polymers do not show that sign change. The first (second) normal stress 
difference for isotropic liquid crystal polymers is positive (negative) and increases (decreases) with increasing shear 
rate. In contrast, active suspensions deep in the isotropic phase can change the sign of the normal stress differences. 
In Fig. [3] the dashed red line shows the dependence of Ni and N 2 on the shear rate for passive liquid-crystal polymers 
(i.e., A = C = 0). 

For contractile suspensions there is a parameter range where the sign of the first normal stress difference changes at 
low shear rates. The upper panel of Fig U shows the relevant parameter range in the A* — £*-plane. The red dashed 
line divides the plane in two regions, one of shear thinning and one of shear thickening. The sign change of Ni appears 
in the shear thinning regime only. The relevant region growth with increasing A* consistent with the limiting case of 
passive nematic liquid crystal polymers (A*— > 1). 

The second normals stress difference for passive suspensions is negative for all shear rates (red dashed line, Fig. 
[3]). However, for active suspension there is a small parameter regime in the shear thickening region Q that shows a 
change in sign as shown in Fig [3J Outside the relevant parameter regions the normal stress difference have the same 
qualitative behaviour as their passive counterparts. The difference is a shift that depends on both parameters A* and 

c*. 

The transition line from shear thickening to shear thinning in Fig. [His not sharp. For a fixed A* value and increasing 
parameter, from = —0.6, a smooth crossover from thickening to thinning in the vicinity of the dashed red line in 
Fig 2] (upper panel) can be observed. The smooth crossover is characterized by a small overshoot that growth and 
disappears as illustrated in the lower panel of Fig. |4j 

In the limit of zero shear, both normal stress differences are vanishing. However, the ratio N2/N1, given by 



hm — 

r^o Ni 



1 



hm — , 

r^o 



is a constant. The so-called viscometric functions are related to the normal stress differences by 



(25) 



(26) 



The ratio of the viscometric coefficients for k ss 0.4 in the passive limit is estimated to ^2/^1 ~ —0.1 which is a 
value typical for polymeric fluids |47j |. The corresponding value for active fluids can vary significantly. For example, 
at C* = 0.3 one obtains * 2 / v I / i ~ 0.2 and for C* = -0.3, ^2/^1 ~ -0.175. This is an additional example for the 
different rheological behaviour of active matter. 



IV. CONCLUSION 



The analytical approach of our simplified model is able to describe some significant rheological features found earlier 
in numerical investigations of more involved models flil Il7l ] . One important feature concerns the effect of activity on 
the zero-shear viscosity which can be significantly larger (contractile) or smaller (extensile) as compared to the passive 
case. Therefore, and since the high-shear viscosity is essentially independent of activity, the intermediate shear-rate 
regime of active suspensions is characterized by either stron g sh ear thinning (contractile) or thickening (extensile). 
Recently, strong shear thinning was experimentally observed [49( in suspensions of Chlamydomonas Reinhardii. 

Furthermore, we investigated the normal stress differences and found a remarkable sign change absent in the passive 
case. This novel effect is characterized by the activity of the suspension. We predict the sign change for the first 
normal stress difference for a wide parameter rage. For example Chlamydomonas Reinhardii may show the predicted 
effect due to strong shear thinning and flow aligning effect similar to rod-like passive suspensions. Moreover, we have 
shown that the corresponding ratio of the viscometric functions deviates markedly from the passive case. 
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FIG. 3: (Color online) Top: The first normal stress difference changes its sign for extensile (£* = 0.6) active particle suspensions. 
The dotted line illuminates negative values of Ni. The curve for contractile suspensions is shifted with respect to the curve of 
passive suspensions (red dashed line). Bottom: The sign change of the second normal stress difference appears for contractile 
suspensions. The dotted line illuminated the positive values of N2. For extensile suspensions there is no sign change, but a 
shift with respect to the passive curve (red dashed line). The remaining parameters of both figures are shown in table [TJ 

In the present study we have not considered the impact of spatially inhomogeneity on the rheology of active matter, 
although this may clearly be relevant [111 ]. Here we have rather assumed that inhomogeneities are small and do not 
change our predictions qualitatively. Indeed, we have performed some preliminary test of our predictions, using an 
one-dimensional, small sine-shaped inhomogeneity. The resulting equations are rather involved and thus not given 
here. An explicit numerical investigation is more reasonable and will be presented elsewhere. Note, that in principle 
the active parameters may also spatially dependent [48j and probably not independent of each other. However, to our 
knowledge there is no rigorous derivation of the parameters from the micro-swimmers that includes the swimming 
mechanism. 

In our study we assumed that the effective interaction between active particles lead to isotropic and nematic order 
as, for example, migrating cells |9J. Similar to the hydrodynamic theory of nematic liquid crystals we coupled the 
generalized Navier-Stokes equations to the second rank symmetric tracclcss order parameter Q. However, in dilute 
suspensions when self-driven particles show swarming motion, the symmetry of the particles oricntational distribution 
leads to a vectorial order parameter P rather than a Q tensor. The rheology for active suspensions related to a 
polarization vector P are studied by Giomi et. al. (l6| . Similar to our findings the model predicts shear thinning, 
thickening and vanishing viscosity. The effect on normal stress differences is not investigated, but we expect no 
significantly different results. In a further extension of the theoretical description a coupled vector-tensor model used 
for passive polarized nano-rod suspensions may show additional rheological effects [HH [52| . 

Very recently, Saintillan used a Fokker-Planck approach to model the orientational dynamics of micro-swimmers in 
extensional flows. He models active stress contributions by force dipoles that are generated by swimming strokes of 
the individual swimmers. Depending on the extensional rate the suspension shows thinning, thickening and a negative 
effective viscosity similar to our results. It is remarkable that, despite different equations used in the passive part of 
the pressure tensor and of the orientational dynamics, all models lead to qualitatively the same rheological behavior. 
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FIG. 4: (Color online)Top: Rheological diagram of active suspensions in the plane spanned by the two activity parameters. 
The red dashed indicated the activity parameter £ c that divides the plane in a shear thinning and a shear thickening region. 
For shear thickening there is a region that shows sign change of the first normal stress difference. On the other hand, in the 
shear thinning regime there exist a parameter range for sign change of the second normal stress difference. The parameter set: 
A = 0.2 and[D 

Bottom: The apparent shear viscosity close to £ c shows a shear thickening and thinning behaviour. The parameter value of 
A = 0.7. 



We expect that the active stress coupling modeled by force dipoles are responsible for the effects discussed here. In 
that sense shear thinning, thickening, zero effective viscosity and sign change in the stress differences are general 
effects in active soft matter. 

We finish with a remark on the relaxation times. The relaxation times of the passive suspension are restricted due 
to Onsager's symmetry relations, viz tq > 0, r p > and tqt p > tq p . On the other hand active particle suspensions 
always dissipate energy due to the process of metabolism. Therefore, a general description needs additional variables 
that are related to the process of metabolism and motility of the swimmers. As a consequence, the relations for the 
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relaxation times in the passive case may be violated. 
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